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Abstract 

Chiral perturbation theory predicts that in quantum chromodymamics light dynamical quarks 
suppress the topological (instanton) susceptibility. We investigate this suppression through direct 
numerical simulation using the Asqtad improved lattice fermion action. This action holds promise 
for carrying out nonperturbative simulations over a range of quark masses for which chiral pertur- 
bation theory is expected to converge. To test the effectiveness of the action in capturing instanton 
physics, we measure the topological susceptibility as a function of quark masses with 2 + 1 dynam- 
ical flavors. Our results, when extrapolated to zero lattice spacing, are consistent with predictions 
of leading order chiral perturbation theory. 

Included in our study is a comparison of three methods for analyzing the topological susceptibil- 
ity: (1) the Boulder hypercubic blocking technique with the Boulder topological charge operator, 
(2) the more traditional Wilson cooling method with the twisted plaquette topological charge op- 
erator and (3) the improved cooling method of de Forcrand, Perez, and Stamatescu and their 
improved topological charge operator. We show in one comparison at nonzero lattice spacing that 
the largest difference between methods (1) and (2) can be attributed to the operator, rather than 
the smoothing method. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.39.Fe 



2 



I. INTRODUCTION 



Chiral perturbation theory predicts the behavior of the topological susceptibility in the 
limit of small quark mass. For improved fermion actions such as Asqtad[l] that lack complete 
chiral symmetry at nonzero lattice spacing, reproducing this prediction is a particularly 
challenging test. 

The Asqtad improvement adds local three-, five-, and seven-link terms to the standard 
staggered fermion action to eliminate tree-level lattice artifacts to order a 2 [1]. Here we 
use it in conjunction with the one-loop improved Symanzik gauge action. This action has 
proven to be highly successful in determining the masses of the light hadrons [2] and a 
variety of quarkonium masses and meson decay parameters [3]. It appears that the most 
satisfactory agreement with experimental values is achieved so far for those quantities that 
are well behaved in the chiral limit. 

Whether improvement is successful in reducing lattice artifacts clearly depends on the 
observable. The Asqtad quark-gluon vertex is not as smooth as that of the more elaborate 
HYP action [4] and zero modes are not treated as rigorously as with the more expensive 
domain wall and overlap actions [5, 6]. Through a rougher vertex, quark propagation might 
be influenced by small instanton-like dislocations. With imprecise zero modes, at small 
quark mass the fermion determinant may fail to suppress adequately configurations with 
nonzero topological charge. 

The gluonic measurement of topological charge by summing the charge density is sensitive 
to the choice of both the discretization of the charge density operator and the smoothing 
or cooling method. Consequently, we found it instructive to compare three methods for 
measuring the charge: 

1. the Boulder definition of the topological charge density [7] with smoothing through 
hypercubic blocking [8] 

2. the more traditional combination of measuring the charge density through the twisted 
plaquette operator and cooling by minimizing the Wilson action [9, 10], and 

3. measuring the charge density with the five- loop improved operator and cooling with a 
five-loop improved gauge action [11]. 
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Throughout, we shall abbreviate these methods with "Boulder/HYP" , "TwPlaq/Wilson" , 
and "5Li/5Li", respectively. 

No previous study of the topological susceptibility has shown satisfactory agreement with 
the predictions of chiral perturbation theory at quark masses much smaller than the strange 
quark mass. Until recently [12, 13, 14], even the expected suppression of the susceptibility 
at small dynamical quark mass has been difficult to detect [15, 16]. We argue that a 
combination of improvements in the lattice action, the smoothing (cooling) technique, the 
topological charge operator, and an 0(a 2 ) extrapolation to the continuum lead to plausible 
agreement with lowest order chiral perturbation theory for small quark masses. 

This article is organized as follows. In Sec. II we introduce notation and review the 
predictions of chiral perturbation theory. Three methods for measuring the topological 
charge are compared and discussed in Sec. III. Results are presented in Sec. IV and discussed 
in Sec. V and conclusions are given in Sec. VI. 

Preliminary results of this study were reported at the Lattice 2002 conference [17]. 



II. TOPOLOGICAL SUSCEPTIBILITY 



The topological charge Q is the integral of the topological charge density p, which is in 
turn defined in terms of the Euclidean color gauge field F° v and its dual JP" , through 



On a lattice of Euclidean space-time volume V the topological susceptibility is the mean 
fluctuation in the topological charge, 

Chiral perturbation theory predicts [18, 19] that at large (Q 2 ) (i.e. large values of mVH, the 
product of the lightest quark mass, the lattice volume, and the chiral condensate parameter), 
and at sufficiently small quark masses, the topological susceptibility is related to the quark 
masses through 

= S 

X (l/m u + l/m d + l/m s y 1 ] 

For three flavors with m u = rrid = m u ^ we may use the PCAC relation between the up and 
down quark mass, pion decay constant, and the chiral condensate to write the susceptibility 



as 

X 4(1 + m u , d /2m s ) ' lJ 
showing that it vanishes linearly in the square of the pion mass in the chiral limit. 

At infinite quark mass (lattice quenched approximation) the susceptibility is finite, requir- 
ing negative curvature corrections in Eq (4), leading asymptotically to a constant [14, 19, 20]. 
Chiral perturbation theory, however, is not expected to converge for masses greater than 
m s , so it provides no guidance there. 

We test Eq (4) against our measured susceptibility, pion mass [21], and pion decay con- 
stant [3] over a range of light quark masses. 

III. MEASURING THE TOPOLOGICAL CHARGE ON THE LATTICE 

In this section we compare three methods for measuring the topological charge. All three 
methods first smooth out ultraviolet fluctuations and then measure the topological charge 
density with a local discretized operator. All are equivalent in the continuum limit. However, 
some perform better at nonzero lattice spacing. To understand the comparison it is first 
useful to briefly review the effects of discretization on the topological susceptibility. 

A. Discretization effects 

The topological susceptibility measured on the lattice differs from the continuum value 
because of discretization errors in the measurement process and in the lattice action itself. 
The principal errors introduced in the measurement process are these: 

1. Instantons with a core size of the order of or less than the lattice spacing are excluded 
altogether, 

2. Small instantons shrink and are erased by prolonged smoothing, 

3. Intermediate sized instantons have a topological charge less than unity, owing to the 
discretization of the charge density operator, 

4. There are 'dislocations' on the lattice: ultraviolet fluctuations that masquerade as 
topological charges, 
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The principal errors introduced by the action are these: 

1. For unquenched gauge configurations the staggered Dirac matrix has inexact zero 
modes so fails to fully "see" and suppress topological fluctuations adequately, and 

2. Lattice artifact flavor (taste) symmetry breaking fails to account for light quark flavors 
correctly. 

B. Comparison of operators and smoothing methods 

We consider the following topological charge density operators: 

1. TwPlaq: The original twisted plaquette operator, defined as an eight-link path with 
a displacement sequence {x, y, —x, —y, z, t, —z, —t} plus rotations. 

2. 5Li: The five-loop improved operator of de Forcrand, Perez, and Stamatescu [11], built 
from a linear combination of five operators in the form of the twisted plaquette, but 
with the plaquettes replaced by various m x n rectangular Wilson loops. 

3. Boulder: The lattice approximation developed for SU(2) by DeGrand, Hasenfratz, 
and Kovacs [7] and refined for SU(3) by Hasenfratz and Nieter [22]. It involves 
a combination of two contorted Wilson loop operators in the fundamental and ad- 
joint representations of SU(3), both defined on closed ten- link paths described by 
unit lattice vector displacements in the sequence {x, y, z, —y, —x, t, x, —t, —x, —z} and 
{x, y, z, —x, t, —z, x, —t, —x, —y} plus rotations and cyclic permutations. This opera- 
tor was optimized to reduce lattice corrections for small instantons with radii close to 
the lattice spacing R « a. 

All three operators are equivalent in the continuum limit, but they are subject to different 
discretization effects. Reference [23] makes a comparison of methods 2 and 3. 

We first investigate how these operators perform with artificial instantons. We created a 
series of gauge configurations containing a single instanton of varying radius and measured 
the topological charge with each operator. Results are plotted in Fig. 1. As expected 
[7] the twisted plaquette operator tends to underestimate the topological charge for small 
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FIG. 1: Topological charge on artificial single instanton configurations as a function of instan- 
ton radius, comparing three observables: the traditional twisted plaquette operator, the Boulder 
operator [7] and the 5Li operatorfll]. 

instantons. The 5Li operator does better. The optimization of the Boulder operator is 
apparent. 

Traditional cooling methods include minimizing the Wilson action [10] and minimizing 
improved actions, such as the 5Li action [11], that provide better scaling as a function of 
instanton radius, so are less likely to erase small instantons. Minimization is done through 
a series of standard heatbath updates (cooling sweeps) at small gauge coupling. Hypercubic 
smoothing [8] was designed to be gentle and local so as to produce a smooth configuration 
with minimal distortion of the topology [24]. The smoothing process involves a series of 
APE blocking steps [25], constrained to lie entirely inside the hypercubes connected to the 
link being smoothed. We use the smoothing coefficients optimized in Ref. [22]. 

To see how cooling or smoothing affects the artificial instantons, we processed them using 
these methods. For the twisted plaquette operator we cooled with ten Wilson gauge action 
updates, for the 5Li operator, ten 5Li updates. For the Boulder operator we smoothed with 
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FIG. 2: Same as Fig. 1, but after smoothing. 

three hypercubic blocking sweeps. The number of smoothing steps in each case was chosen 
for stability of x under further smoothing, as we shall discuss below. Results are shown in 
Fig. 2. It is evident that small instantons preserve their topological charge best with the 
Boulder /HYP method. 

Next we examine how the three methods perform on one of the gauge ensembles in our 
study, namely the 20 3 x 64 a = 0.12 fm set with quark masses am u ^ = 0.01 and am s = 0.05 
[27]. We measured the topological susceptibility as a function of cooling or smoothing step 
and compared the results in Fig. 3. We see that it is reasonable to read off the Boulder/HYP 
susceptibility after three HYP sweeps and the TwPlaq/ Wilson and 5Li/5Li susceptibilities 
after ten cooling sweeps. We made these arbitrary choices in an effort to compromise between 
preserving small instantons and reaching stability in the observable. We have tested these 
choices in a few cases and find within statistical errors that our results are insensitive to 
increasing these values by a factor of two or three. See also Ref. [23]. 

It is clear that at this lattice spacing the TwPlaq/Wilson method gives a lower suscep- 
tibility than the other methods. The 5Li/5Li result is closer to but still lower than the 
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FIG. 3: Topological susceptibility as a function of smoothing or cooling step on the 20 3 x 64 
dataset with am u ^ = 0.01 and am s = 0.05, comparing four techniques: the Boulder topological 
charge operator with HYP smoothing [8], the twisted plaquette operator [9] with both Wilson 
action cooling [10] and HYP smoothing, and the 5Li operator with 5Li cooling [11]. (Results 
are expressed in units of ro, the Sommer parameter [26]). The arbitrary scale conversion counts 
three cooling steps for one HYP smoothing sweep. Susceptibilities are measured on subvolumes as 
explained in Sec. IV. 

Boulder/HYP result. To determine to what extent the difference in TwPlaq/Wilson is at- 
tributable to the observable and what to the cooling method, we also measured the TwPlaq 
susceptibility on HYP smoothed lattices. The result (TwPlaq/HYP) shown in Fig. 3 is quite 
close to the TwPlaq/Wilson result. So the choice of operator appears to account for the 
largest discrepancy. 
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We have also measured the topological susceptibility on the companion quenched 20 3 x 64 
a = 0.12 fm ensemble using three methods and found similar inequalities: TwPlaq/Wilson 
gave x4 = 0.036(2); 5Li/5Li, 0.051(3); and Boulder/HYP, 0.054(3). 
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FIG. 4: Topological susceptibility vs pion mass squared in units of ro on the a = 0.12 fm lattices, 
comparing the twisted plaquette plus Wilson method method and the Boulder plus hypercubic 
blocking method. Also shown are two results for the 5Li/5Li method. The solid line shows the 
prediction of leading order chiral perturbation theory. The quenched result is shown at the extreme 
right. The susceptibility is measured on subvolumes, as explained in Sec. IV. 

Finally we considered the entire 20 3 x 64 dataset with a ~ 0.12 fm tabulated in Ta- 
ble I. Taking ten cools and three HYP sweeps for the comparison, we plot the result in 
Fig. 4. Throughout the entire mass range the TwPlaq/Wilson susceptibility is about 2/3 
the Boulder/HYP value. 
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Since all three methods are expected to give the same continuum limit, as they do for 
quenched QCD with the Wilson plaquette action [28], the discrepancy we observe at lattice 
spacing a = 0.12 fm must be due lattice artifacts. To decide what the continuum limit is 
and which method is closer to it, one should do a detailed scaling study. We present results 
of a partial study in Sec. IV, but here attempt to interpret the differences based on our 
observations on smooth artificial instantons: 

These discrepancies are of a magnitude that would be expected from our measurements 
of the charge of artificial instantons of intermediate size. The average instanton radius is 
expected to be approximately 0.3 fm [29] or 2.3a on these a = 0.12 fm lattices with significant 
contributions from radii as small as one lattice unit. From Fig. 2 where ideally Q = 1 we 
see that the TwPlaq/Wilson method underestimates the charge by Q(2.3a) = 0.85 after 
cooling, so one would expect an underestimate of about 0.85 2 = 0.7 in the susceptibility for 
the average instanton. For smaller instantons the TwPlaq/Wilson method gives Q 2 (1.5a) = 
0.03 after cooling, instead of one. From Fig. 3, we find a ratio of 0.58(4) between the 
TwPlaq/WilsonlO and Boulder/HYP3 susceptibilities at (m^ro) 2 = 0.738. The ratio is 
approximately the same throughout the entire mass range. For the 5Li method we have 
Q 2 (2.3a) = 0.96 and Q 2 (1.5a) = 0.80. By comparison the ratio 5Li/5Li to Boulder/HYP3 is 
0.87(7) at (m n ro) 2 = 0.738. Consequently, one may wonder whether the apparent agreement 
at a = 0.12 fm between the TwPlaq/Wilson method and chiral perturbation theory at 
2 < (m^ro) 2 < 3 is the result of compensating errors. 

IV. RESULTS 

We measured the topological susceptibility using the the Boulder/HYP method on two 
sets of gauge configurations generated with three flavors of light Asqtad quarks of varying 
masses, one set with lattice spacing approximately 0.12 fm throughout and the other, 0.09 
fm [27]. The corresponding matched quenched configurations are also included. The data 
sample is tabulated in Table I. 

Besides measuring the susceptibility on the entire lattice volume, we increased our statis- 
tics by measuring on smaller subvolumes. The probability distribution follows a Gaussian 
in Q with width proportional to the volume. The width is decreased as the volume is 
decreased, leading to the same relative error in the determination of the susceptibility for 
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the same sample size on the smaller volume. To the extent that the subvolume measure- 
ments are uncorrelated, the sample is effectively increased by a factor equal to the number 
of subdivisions, so the error should decrease by the square root of this number. This pro- 
cess cannot be continued indefinitely, since we should eventually discover strong correlations 
among adjacent subvolumes. 

Measuring the susceptibility on subvolumes may even be indicated for actions and up- 
dating algorithms that give persistent global charge, but fluctuating local charge densities. 

Accordingly, we divided the lattices along the time dimension into three hypercubic sub- 
volumes mostly separated by a small unused space. For the 28 3 x 96 lattices these subvolumes 
of size 28 4 were constructed from imaginary time ranges [0,27], [32,59], and [64,91], and 
for the 20 3 x 64 lattices, of size 20 4 from time ranges [0,19], [22,41], and [44,63]. Since 
the boundary condition on the subvolume is not periodic, this practice splits instantons. 
However, the rms charge in the subvolumes (in the range 4 < \Q\ < 7) seems large enough 
to ensure that the boundary effects are not significant. 

We measured the correlations in the topological charge history between charges in differ- 
ent subvolumes. That is, if we denote by Q u k the charge measured in subvolume u on gauge 
configuration A; in a data sample with N gauge field configurations, we define the correlation 
coefficient to be 

Cuv = J7^2QukQvk/(\Qu\\Qv\) (5) 
k 

where \Q U \ = J (Qt) * s the rms charge on the subvolume u. The 33 correlation coefficients 
for our entire data set are roughly Gaussian distributed about zero with a mean of —0.004 
and width of 0.09. Thus we feel confident that we may treat the subvolume measurements 
as statistically independent observations. 

The measurements are clearly correlated in Monte Carlo time. The Asqtad dynamical 
lattices were saved every sixth molecular dynamics trajectory. The 20 3 x 64 quenched lat- 
tices were saved every tenth quasi-heatbath sweep and the 28 3 x 96 every fiftieth. We made 
charge measurements on all available lattices. A particularly striking example is given by 
the time history for the total charge on the 28 3 x 96 lattice with three degenerate quark 
masses am u ^ = am s = 0.31, as shown in the upper panel of Fig. 5 [17]. The horizontal scale 
counts molecular dynamics trajectories. Time histories for the subvolume charges for the 
same dataset are shown in Fig. 6 where the prominent oscillations are much less evident. 
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Other time histories shown in the lower two panels of Fig. 5 do not show such a striking 
effect. On the 28 3 x 96 datasets the autocorrelation length, measured on the subvolumes 
by summing the autocorrelation coefficient, ranges from approximately 10 trajectories for 
o^u.d — 0.0062 to 35 trajectories at am u ^ = 0.031. Charge measurements on the corre- 
sponding quenched lattices are only slightly correlated. On the coarser 20 3 x 64 lattices 
we find weaker autocorrelation, but still roughly monotonically increasing with quark mass 
from fewer than six trajectories at am u ^ = 0.007 to ten trajectories at am u ^ = 0.050. This 
trend appears to be contrary to some expectations [30]. 

The autocorrelations in topological charge found with the Asqtad action and 2 + 1 flavors 
of quarks appear to be longer than those found with the conventional thin link staggered 
fermion action and two flavors. Shown in Fig. 7 is a comparison of the topological charge 
history from the ensemble of Fig. 5b and an ensemble generated with the conventional 
unimproved thin-link staggered fermion algorithm [31]. These simulations were done at 
approximately the same value of (m n ro) 2 (unimproved 1.06, improved 1.23) and lattice 
spacing (unimproved 0.10 fm, improved 0.09 fm). It is apparent that the configurations 
decorrelate less rapidly with the improved action and extra flavor. 

Statistical errors in the topological susceptibility are determined by taking the larger of 
the error corrected for autocorrelations and the error obtained by extrapolating to infinite 
bin size. 

We have measurements of the topological susceptibility at two lattice spacings. Thus we 
may venture an extrapolation to the continuum limit. This is done by first interpolating 
in (rom^) 2 the topological susceptibility on the coarse lattice to the three pion mass values 
where we have measurements on the finer lattice. We then do a linear extrapolation at fixed 
(m^ro) 2 to zero a 2 jr\ as shown in Fig. 8. We see that the quenched susceptibility rises with 
decreasing lattice spacing, but falls for the unquenched lattices with a slope that increases 
with quark mass. 

We have done a comparable extrapolation using the TwPlaq/Wilson method, but with 
only one of the a = 0.09 fm dynamical ensembles. Results are also shown in Fig. 8. Within 
the limitations of these few points, we find satisfactory agreement between the extrapolated 
values found with both methods. 

Our main results for the topological susceptibility are summarized in Table I and Fig. 9. 
We see that within errors the continuum extrapolation gives reasonable agreement with the 
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FIG. 5: Full volume topological charge after three HYP sweeps vs molecular dynamics trajectory 
for the 28 3 x 96 data set with (a) am Ui d = am s = 0.031, (b) am U: d = 0.0124, am s = 0.031, and (c) 
am U) d = 0.0062, am s = 0.031. In the last case two separate time series are plotted with the second 
starting after the break at 700 trajectories. 

prediction of leading order chiral perturbation theory, Eq (4). 



V. DISCUSSION 

The combined effect of the discretization artifacts discussed in subsection III A is that, 
depending on the method of measurement, the lattice topological susceptibility x is both 
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FIG. 6: Same as Fig. 5a, but measured on three 28 4 subvolumes defined by time slices (a) [0,27], 
(b) [32,59] and (c) [64,91]. 



additively and multiplicatively renormalized relative to the continuum value, \ [9, 33]: 

X(a, m q ) = M(a, m q f x(m q ) + A(a, m q ), 



(6) 



where we understand the susceptibilities to be expressed in some appropriate choice of 
physical units. Having some indication of the scaling behavior of these methods, we suggest 
a scenario for the behavior of the functions M and A in Eq. (6). 

To suppress the ultraviolet fluctuations prior to measuring the topological charge, it is 
necessary to smooth the gauge configurations. In the quenched theory, smoothing eventually 
drives A to 0. However, since fluctuations involving instantons of sizes at or below the cutoff 
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trajectory 

FIG. 7: (a) Topological charge after three HYP sweeps vs molecular dynamics trajectory for the 
conventional unimproved staggered fermion action on a 24 3 x 64 data set with two degenerate 
quark flavors am u d = 0.01, compared with (b) the result from the Asqtad action from Fig. 5(b). 

are always missed, one expects M < 1, at least for operators that do not overestimate the 
charge of small instantons. Therefore, x < X ■ Indeed, we see from Fig. 4 and the 
continuum extrapolation in Fig. 8 that at a = 0.12 fm all three methods underestimate the 
quenched susceptibility. 

In full QCD the virtual fermions screen the topological charge, but they screen everything 
they see, including dislocations. Smoothing removes dislocations and small instantons; those 
that are not involved in the screening of more extended topological charges do not cause 
trouble, but those that are leave behind the more extended (and now less screened) charges 
that are relatively stable under smoothing. This effect increases the susceptibility and shows 
up as an anomalous A that does not vanish with continued smoothing. 

As we take the continuum limit in full QCD, we expect M — > 1, as in the quenched 
case, and also A — > 0, as the lattice action, combined with smoothing, gradually suppresses 
dislocations. The two trends act in opposite directions; to decide which dominates, we must 
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FIG. 8: Continuum extrapolation of the topological susceptibility found with the Boulder/HYP 
and for the lightest quark mass the TwPlaq/ Wilson methods. 

consider the quark mass dependence. There are two regimes. The first is for large m q , where 
M dominates and A is small. In this "instanton-dominated" region, the behavior is similar 
to the quenched theory, and discretization effects lead to an underestimate x(a, m q ) < x( m g) 
at nonzero a. On the other hand, at small quark mass the continuum susceptibility is small. 
The factor M is expected to depend only weakly on m q , so the lattice measurement is dom- 
inated by A. In this "dislocation-dominated" regime there is likely to be an enhancement 
x(a,m q ) > x{ m q) & t nonzero a. Indeed, from Fig. 4 and the continuum extrapolations at 
(m^ro) 2 = 0.646 and infinity in Fig. 8 we see that at a = 0.12 fm all three methods overesti- 
mate the susceptibility expected at the nearby point (m^ro) 2 = 0.738 and underestimate it 
at infinite quark mass. With the Boulder/HYP method, our results suggest that the scaling 
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FIG. 9: Topological susceptibility vs pion mass squared in units of ro with dynamical Asqtad 
quarks. The solid line shows the prediction of leading order chiral perturbation theory. The dotted 
line gives a phenomenological proposal for nonperturbative behavior [32]. The quenched result is 
shown at the extreme right. 

slope decreases as the quark mass is increased. 

Given that the mean gauge action decreases by orders of magnitude under smoothing 
at a ~ 0.1 fm, it is reasonable to assume that dislocations are abundant and participate 
significantly in screening more extended topological charges. Erasing them then contributes 
to A. The number of instantons is not a strong function of quark mass. (The Wilson gauge 
action after twenty cooling sweeps is dominated by instantons, so gives a measure of their 
number. We find that the mean action density on the coarse lattices varies by less than 10% 
across the dynamical and quenched ensembles in this study). We thus expect A to vary 
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at most weakly with the quark mass in the dislocation-dominated region of quark mass. 
Our results for the a = 0.12 fm lattices plotted in Fig. 4 show that for all methods the 
susceptibility levels off below (m n r Q ) 2 ~ 2. 

But leveling off can also be attributed to shortcomings of the lattice fermion formulation 
itself. The would-be zero modes of the staggered lattice Dirac matrix are not at exactly zero. 
If the deviation is comparable to the virtual quark mass, some of the topological modes are 
not properly screened. Consequently, x(a, m q ) does not vanish as it should as the quark 
mass is reduced. The chiral limit is then governed by the dependence of A on m q . 

VI. CONCLUSIONS 

We have measured the topological susceptibility on lattices generated with Asqtad im- 
proved staggered fermions of varying mass and two lattice spacings, a = 0.12 fm and a = 0.09 
fm. We have compared three methods for measuring the topological charge on these lattices 
and selected the Boulder method with hypercubic blocking, since it appears best capable at 
these lattice spacings of preserving small instantons. We show in one comparison at a = 0.12 
fm that the largest difference between the Boulder/HYP and TwPlaq/Wilson methods can 
be attributed to the operator, rather than the smoothing method. We find that at both 
lattice spacings there is clear evidence that dynamical quarks suppress topological fluctua- 
tions and the suppression increases with decreasing quark mass and with decreasing lattice 
spacing. However, at fixed nonzero lattice spacing, lattice artifacts at small quark masses 
are still substantial, and the susceptibility does not decrease as expected from chiral per- 
turbation theory at small quark masses. Nevertheless, an 0(a 2 ) extrapolation of our data 
to zero a gives results that are consistent with the leading order prediction. Within the 
limitations of our statistics, this is the first study to show satisfactory agreement with the 
predictions of chiral perturbation theory at quark masses much smaller than the strange 
quark mass. 
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TABLE I: Topological susceptibility vs pion mass squared in units of tq. The result xv is t ne 
susceptibility computed on the full lattice volume. The result x 1S computed on three subvolumes 
for improved statistics. 
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